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I.  Introduction 


The  Modified  Point  Mass  Trajectory  Model  is  the  primary  method  of  trajectory 
simulation  used  in  the  preparation  of  Firing  Tables.  This  model  requires  three  types 
of  input  data:  projectile  mass  properties,  aerodynamic  coefficients  and  the  performance 
parameters  determined  from  range  testing.  This  report  discusses  an  initial  attempt  to 
determine,  for  trajectory  modeling,  the  aerodynamic  drag  of  the  155mm,  DPICM,  M864 
base-burn  projectile  (Figure  1)  from  the  HAWK  Doppler  radar  data.^^  A  detailed  treatise 
of  base-burn  projectile  ballistic  modeling  technology  is  reported  by  Niles-Erik  Gunners, 
Kurt  Andersson  and  Rune  Hellgren  in  Reference  5,  Chapter  16,  “Base-Bleed  Systems  for 
Gun  Projectiles”. 


II.  Reduction  of  HAWK  Doppler  Radar  Velocimeter  Flight 

Data 

The  frame  of  reference  for  all  vectors  is  a_ground-fixed,  orthonormal,  right-handed 
Cartesian  coordinate  system  with  unit  vectors  ( 1,  2  and  3 ).  The  1  axis  is  the  intersection 
of  th^  vertical  plane  of  fire  and  the  horizontal  plane  and  pointing  in  the  direction  of  fire. 
The  2  axis  is  parallel  to  the  gravity  vector,  g,  and  opposite  in  direction.  The  3  axis 
completes  the  right-handed  coordinate  system. 

The  slant  range  rate  of  change  as  measured  by  HAWK  Doppler  radar  is  recorded  on 
a  digital  tape.  The  first  step  is  to  smooth  the  data  and  determine  the  time  derivative. 
Least  squares  fits  (second-degree  polynomials)  to  the  data  are  determined  for  0.56  second 
intervals  (fifteen  point  smoothing)  along  the  trajectory.  The  slant  range  rate  of  change  (r) 
and  time  derivative  of  the  slant  range  rate  of  change  (f )  are  obtained  from  the  quadratic 
fit  at  the  midpoint  of  the  fifteen  point  interval. 

An  estimated  trajectory  is  generated  separately  using  the  projectile  mass  properties, 
launch  data,  atmospheric  conditions,  estimated  aerodynamic  coefficients  and  estimated 
base  drag  reduction  factor  during  base-burn  motor  burning  as  a  function  of  time.  The 
trajectory  is  adjusted,  using  factors  on  base  drag  reduction  during  motor  functioning  and 
lift,  to  match  the  observed  impact  data.  A  trajectory  velocity  (tt^. )  is  calculated  using  the 
HAWK  radar  smoothed  slant  range  rate  of  change  (r)  and  the  estimated  trajectory  slant 
range  rate  of  change  (r't )  and  velocity  (tTt )  as  follows: 

Ur  =  (r/rt)tr,  (1) 


where: 


rf  =  Ut  cos  (r, ,  u,)  =  (  rj  •  tT, )  /  r, 
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Note:  Subscript  i  refers  to  quantities  determined  from  the  estimated  trajectory  and  those 
with  subscript  r  are  obtained  using  both  the  HAWK  radar  data  and  the  estimated 
trajectory. 

A  trajectory  acceleration  (tfl.)  was  calculated  using  the  time  derivative  of  the  HAWK  radar 
slant  range  rate  of  change  (f)  and  the  estimated  trajectory,  time  derivative  of  the  slamt 
range  rate  of  change  (ft)  and  acceleration  (tit ),  using  the  following  two  formulations: 

Ur  =  (r  /rt)  ut  +  [(r’f  f  -  r  r't)  /  f(*]  u,  (2) 

and 

Ur  =  (f/ft)Ut  (3) 


where: 


-  (n  *  )  rt  } /rt^ 

and 


The  mean  of  the  results  were  similar  for  both  of  the  Ur  representations,  however,  the 
variation  (spread)  of  the  results  were  significantly  improved  using  equation  3  and  it  was 
used  for  determining  the  results  presented. 


III.  Determination  of  Aerodynamic  Drag  Coefficient 

The  mass  of  the  projectile,  atmospheric  conditions,  estimated  trajectory  data  and  the 
Doppler  slant  range  rate  of  change  and  time  derivative  of  the  slant  range  rate  of  ch^Lnge 
provide  the  necessary  inputs  to  determine  the  aerodynamic  drag.  The  following  inverse 
solution  of  the  point-mass  equations  of  motion  is  then  used  to  compute  the  aerodynamic 
drag  (Cor)- 


Cpr  =  -  ((tTr  -«?)•(«;  -  g  -  A)]  8  m/(7r  p  d®  V®  )  (4) 
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IV.  Determination  of  Base  Drag  Reduction  Factor 


The  base-burn  motor  is  intended  to  increase  the  range  of  the  projectile  by  reducing 
the  base  drag.  An  initial  estimate  of  the  aerodynamic  zero  yaw  drag  force  coefficient  {Cdo  ) 
and  the  base  drag  component  {Cdb  )  was  determined  for  the  M864  base-burn  projectile, 
with  an  inert  base-burn  motor,  using  the  semi-empirical  drag  estimation  model  known  as 
McDrag.®  A  base  drag  reduction  factor  (/bzj)  was  then  defined  to  quantify  the  amount  of 
base  drag  reduction  as  follows: 


fBD  =  (C/Jo  -  CbJ/Cbb,  (5) 

A  base  drag  reduction  factor  of  unity  indicates  that  all  the  estimated  base  drag  is  ehmi- 
nated,  whereas  a  value  of  zero  means  none  of  the  estimated  base  drag  is  eliminated.  A  fsD 
greater  than  one  indicates  that  thrusting  is  present,  since  more  than  just  the  base  drag 
has  been  eliminated. 


V.  Results 

The  aerodynamic  drag  results,  obtained  from  the  procedure  described  above,  are 
presented  for  a  sampling  of  the  M864  projectiles,  including  three  M864  projectiles  with 
inert  base-burn  motors.  These  projectiles  were  fired  at  Yuma  Proving  Ground,  AZ  during 
May  1987.  The  base-burn  motors  were  designed  to  burn  out  at  approximately  30  seconds. 
The  aerodynamic  drag  coefficients  for  the  projectiles  with  inert  motors  are  presented  as 
a  function  of  Mach  number  and  for  the  projectiles  with  live  motors  base  drag  reduction 
factors  are  presented  as  a  function  of  time  of  flight. 

The  sampling  includes  projectiles  with  inert  and  live  base-burn  motors  fired  with 
propelling  charges:  M4A2,  charge  7W;  M119A2,  charge  7R;  and  M203E2,  charge  8R.  The 
sampling  contains  an  inert  base-burn  motor  fired  at  a  quadrant  elevation  of  748  mils  and 
live  base-burn  motors  fired  at  three  quadrant  elevations  (approximately  500,  750  and  1150 
mils)  with  each  of  the  propeUing  charges. 

Figures  2  through  4  present  the  aerodynamic  drag  coefficients  versus  Mach  number 
for  projectiles  with  inert  base- burn  motors  fired  with  each  of  the  propelling  charges  at 
a  quadrant  elevation  of  748  mils.  Figures  5  through  13  present  the  base  drag  reduction 
factors  versus  time  of  flight  for  projectiles  with  hve  base-burn  motors  fired  with  each  of 
the  propelling  charges  at  quadrant  elevations  of  approximately  500,  750  and  1150  mils. 

Figures  2,  3  and  4  show  excellent  agreement  of  the  aerodynamic  drag  coefficients 
for  the  inert  base-burn  motor  projectiles  fired  with  the  three  propelling  charges.  Fig¬ 
ures  5  through  13  show  a  good  correlation  of  the  base  drag  reduction  factor  with  time 
of  flight  for  the  hve  base-burn  motor  projectiles.  The  base  drag  is  consistently  reduced 
by  approximately  50  percent  and  the  results  show  an  increase  in  this  percentage  for  the 
higher  quadrant  elevations  indicating  a  possible  correlation  with  local  atmospheric  pres- 
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sure.  There  is  some  irregularity  in  the  base  drag  reduction  factor  for  the  transonic  velocities 
(Mach  numbers:  .95  to  1.05).  This  is  identified  where  evident  on  some  of  the  figures.  The 
irregularity  is  probably  due  to  the  error  in  the  transonic  aerodynamic  inputs  and/or  the 
Mach  number  determined  from  the  HAWK  radar  data.  The  base  drag  reduction  factors 
also  indicate  a  delay  in  base  drag  reduction  at  the  beginning  of  approximately  one  second 
and  a  base-burn  motor  burnout  time  of  approximately  33  seconds. 

The  large  variation  in  the  base  drag  reduction  factor  for  the  high  quadrant  elevation 
rounds  is  probably  due  to  a  combination  of  the  following:  the  magnitude  of  the  total 
drag  coefficient,  the  data  reduction  methodology  including  the  estimated  trajectory,  the 
estimated  aerodynamic  coefficients  and  the  HAWK  Doppler  radar  capability. 


VI.  Analysis  of  Aerodynamic  Drag  Results 


An  aerodynamic  drag  coefficient  was  determined  for  the  M864  projectile  based 
on  the  inert  base-burn  motor  projectile  results.  In  addition,  an  estimated  inert,  base  drag 
component  was  determined  for  the  M864  projectile  with  an  inert  base- burn  motor 

by  subtracting  the  McDrag  estimates  of  head,  skin  friction,  band  and  boat  tail  drag  from 
the  determined  Con  for  the  inert  M864  projectile.  Figure  14  presents  the  determined  Cd„ 
and  Cdb,  functions  of  Mach  number. 

The  base  drag  reduction  factors  for  times  of  flight  between  three  and  33  seconds  were 
analyzed  using  a  stepwise  multiple  regression  technique. r  The  base  drag  reduction  factors 
were  fit  as  a  function  of  time  of  flight  and  the  difference  in  a  reference  (standard)  air 
pressure  (Pr)  ,  1013.25  mb,  and  the  local  air  pressure  (P)  divided  by  P,  to  account  for 
the  apparent  variation  with  quadrant  elevation.  A  review  of  the  residuals  as  a  function 
of  Mach  number  suggested  an  additional  function  of  Mach  number,  fi  ,  which  forms  a 
modified  base  drag  coefficient,  {Cpg.  /i),  for  the  M864  base-burn  projectile.  Figure  14 
presents  the  modified  base  drag  coefficient,  (Cdb  fi),  while  Figure  15  presents  the  fi 
factor  as  a  function  of  Mach  number.  The  value  of  {Cdb  /i)  25  percent  higher  for  Mach 

numbers  less  than  0.9  and  11  percent  lower  at  Mach  number  2.0  than  the  theoretically 
determined  Cdb  •  A.  theoretical  discussion  of  drag  reduction  for  base-burn  projectiles  is 
presented  in  Reference  5,  Chapter  16. 

The  stepwise  multiple  regression  fitting  function  including  /i  is  as  follows: 


[{Cd.  -  Cd^)/(Cdb,  fi)]  =  /2  +  /4  (Pr  -  P)/Pr  (6) 

Where: 

•  Factor  /2  is  a  function  of  time  of  flight. 

•  Factor  /<  is  a  constant. 

Note:  A  factor  /j  multiplying  /i  is  added  in  the  final  formulation. 
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The  factor,  /2  ,  as  a  function  of  time  of  flight,  determined  by  the  stepwise  multiple 
regression  process  for  the  nine  M864  base-burn  test  projectiles  (Figures  5  through  13),  is 
presented  in  Figure  16  and  the  constant,  ,  determined  simultaneously,  was  0.30.  The 
residuals  of  the  stepwise  multiple  regression  fitting  process  are  presented  in  Figures  17,  18 
and  19  versus  time  of  flight,  Mach  number  and  local  air  pressure  respectively.  The  root 
mean  square  error  of  residuals  is  0.12  and  the  residuals  show  that  the  fitting  process  has 
removed  the  symmetric  biases  due  to  time  of  flight,  Mach  number  and  locsd  air  pressure. 

Summarizing: 

•  An  aerodynamic  drag  coefficient  (Cdo)  was  determined  for  the  M864  projectile 
based  on  the  inert  base-burn  motor  projectile  results. 

•  A  base  drag  component  {Cdb  )  was  theoreticadly  determined  for  the  M864  pro¬ 
jectile  with  an  inert  base-burn  motor  based  on  minus  the  fore  body  drag 
estimated  by  McDrag. 

•  A  factor,  /i  ,  multiplied  by  the  inert  projectile  base  drag  component  {Cdb  ) 
used  to  represent  an  effective  inert  base  drag  coefficient  for  the  M864  base-burn 
projectile  as  a  function  of  Mach  number. 

•  A  factor,  fi  .  is  used  to  represent  the  effect  of  time  on  drag  reduction  during  the 
burning  phase  of  the  M864  base- burn  projectile. 

•  A  change  in  drag  reduction  with  local  air  pressure  was  assumed  to  account  for 
the  observed  variation  of  the  M864  base-burn  projectile  base  drag  with  quadrant 
elevation.  This  apparent  effect  with  local  air  pressure  was  further  assumed  to 
take  the  form  fii  Pr  —  P  )  /  Pr ,  where  the  factor,  /4  ,  is  a  constant. 

VII.  Exterior  Ballistic  Trajectory  Simulation  Model 

The  reduction  in  base  drag  during  base-burn  motor  functioning  is  shown  to  be  a 
function  of  time  of  flight  and  to  a  lesser  degree  a  function  of  local  air  pressure.  A  constant, 
/a ,  multiphed  by  /j  is  added  to  create  a  simulation  methodology  containing  flexibility  for 
matching  experimental  radar  and  impact  or  air  burst  data.  The  constants,  /a  and  , 
can  then  be  used  effectively  as  parameters  for  matching  the  measured  impact  or  air  burst 
range  firing  data. 

The  combined  base  drag  reduction  term,  —  (  Cog.  /i  )  [  /a  /j  +  f*  {Pr  —  P)  /  Pr],  is 
then  added  to  the  projectile’s  inert  total  drag  force  coefficient  term,  Crh>  +  ((?ae)’, 

of  the  Modified  Point  Mass  Trajectory  Model  for  Rocket  Assisted  Projectiles. 

The  total  drag  computed  by  the  Modified  Point  Mass  Trajectory  Model  for  Rocket 
Assisted  Projectiles  is  then  as  follows: 
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D  = 


fo 

8  m 


{  Co.  -  (  Cd.  h)lhh  +  U(P,  -  P)/Pr]  +  Cn,  {Qa.f  }  V  (7) 


Where: 

•  The  factor,  /o ,  can  be  used  as  a  function  of  quadrant  elevation,  instead  of  the 
ballistic  coefficient,  for  matching  experimented  impact  or  air  burst  range  firing 
data  for  projectiles  without  base-burn  motors.  A  value  of  1.0  was  used  for  fo 
because  this  parameter  was  not  used  for  matching  the  experimental  firing  data  of 
the  M864  ba®e-burn  projectile. 

Note:  The  ballistic  coefficient,  m,  /  ( /o  d*  ) ,  as  a  function  of  quadrant  elevation 
is  currently  in  use  as  the  parameter  for  matching  experimental  impact  or 
air  burst  range  firing  data.  The  use  of  fo ,  previously  identified  as  i ,  has 
the  advantage  of  being  a  nondimensional  factor. 

•  The  term,  (  Cpa  /i  )  (  /2  /a  ) ,  is  used  to  represent  the  drag  reduction  as  a  function 
of  Mach  number  and  time  of  flight  and  for  matching  experimental  impact  or  air 
burst  range  firing  data  of  individual  or  groups  of  projectiles. 

•  The  factor,  fi  ,  as  illustrated  in  Figure  15,  is  used  to  represent  the  change  in 
drag  reduction  as  a  function  of  Mach  number.  /]  was  determined  from  the 
analysis  of  the  HAWK  radar  data. 

•  The  factor,  /j  ,  as  illustrated  in  Figure  16,  is  used  to  represent  the  change  in 
drag  reduction  as  a  function  of  time  of  flight.  /2  was  also  determined  from 
the  analysis  of  the  HAWK  radar  data. 

•  The  factor,  /a  can  be  used  as  a  parameter  for  matching  experimental  impact 
or  air  burst  range  firing  data. 

•  The  term,  (  Cdb  fi)[f4  {Pr  —  P)  /  Pr],is  used  to  represent  the  drag  reduction 
due  to  an  apparent  effect  of  the  local  air  pressure. 

•  The  factor,  /4  ,  is  a  constant  during  base-burn  motor  functioning  and  is  zero 
thereafter,  was  determined  to  be  0.30  based  on  the  analysis  of  the  HAWK 
radar  data. 

•  The  factor,  Q ,  is  a  constant  used  to  compensate  in  part  for  the  approximations 
in  the  Modified  Point  Mass  Trajectory  Model.  A  value  1.2,  normally  used  for 
artillery  projectiles,  was  used  for  this  parameter. 

In  the  above  model  the  nondimensional  factor,  fo ,  replaces  the  ballistic  coefficient 
as  the  parameter  used  for  matching  experimental  impact  or  air  burst  range  firing  data  of 
projectiles  without  base-burn.  In  that  role,  fo  would  normally  be  a  function  of  quadrant 
elevation.  The  factors,  /i  through  /<  ,  should  be  effective  parameters  for  representing  the 
HAWK  radar  results  and  matching  the  experimental  impact  or  air  burst  range  firing  data 
of  the  155mm,  M864  base-burn  projectile. 
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VIII.  Conclusions 


The  aerodynamic  drag  coefficients  determined  for  the  M864  projectiles  with  inert  base- 
burn  motors  show  excellent  agreement  for  the  flight  Mach  numbers  of  projectiles  fired  with 
propelhng  charges:  M4A2,  charge  7W;  M119A2,  charge  7R;  and  M203E2,  charge  8R. 

The  reduction  in  base  drag  during  base-burn  motor  functioning  correlated  very  well 
with  time  of  flight  and  to  a  lesser  degree  with  local  air  pressure.  These  results  support  the 
simplistic  addition  of  a  base  drag  reduction  term  to  the  acceleration  due  to  drag  equation 
of  the  Modified  Point  Mass  Trajectory  Model  for  Rocket- Assisted  Projectiles  for  simulating 
the  ballistic  trajectory  of  the  155mm,  M864  base-burn  projectile. 
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155mm,  DPICM,  M864  Base-Burn  Projectile 

Sketch 


Length  of  Projectile 

calibers 

5.79 

Nose  Length 

calibers 

3.42 

Cylinder  Length 

calibers 

1.86 

Boattail  Length 

calibers 

.50 

Boattail  Angle 

degrees 

3.00 

Mass  Properties 

Mass 

kgs 

46.95 

(lbs) 

103.5 

Mass  of  Fuel 

kgs 

1.21 

(lbs) 

2.67 

Center  of  Gravity 

cm  from  nose 

58.8 

(in  from  nose) 

23.16 

Moments  of  Inertia 

Axial 

kg-m* 

.158 

(Ib-ft*) 

3.75 

Transverse 

kg-m* 

1.657 

(Ib-ft*) 

39.32 

Figure  1.  Physical  characteristics  of  the  155mm,  DPICM,  M864  base-bum  projectile. 
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Aerodynamic  drag  coefficient  versus  Mach  mimber  for  round  number  4376 
dth  inert  base- burn  motor,  fired  with  propelling  charge  M119A2,  7R,  at  a 
luadrant  elevation  of  748  mils. 
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Figure  4.  Aerodynamic  drag  coefficient  versus  Mach  number  for  round  number  4341 
with  inert  base-burn  motor,  fired  with  propelling  charge  M203E2,  8R,  at  e 
quadrant  elevation  of  748  mils. 
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Figure  6.  Base  drag  reduction  factor  versus  time  of  flight  for  round  number  1034  fired 
with  propelling  charge  M4A2,  7W,  at  a  quadrant  elevation  of  750  mils. 
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0.0  5.0  10.0  15.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0 

Time  of  Flight  (Seconds) 

Figure  7.  Base  drag  reduction  facto’*  versus  time  of  flight  for  round  number  1013  fired 
with  propelling  charge  M4A2,  7W,  at  a  quadrant  elevation  of  1150  mils. 
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;  reduction  factor  versus  time  of  flight  for  round  number  1050  fired 
lelling  charge  M119A2,  7R,  at  a  quadrant  elevation  of  750  mils. 
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;  reduction  factor  versus  time  of  flight  for  round  number  4202  fired 
lelline  charee  M119A2,  7R,  at  a  quadrant  elevation  of  1150  mils. 
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Figure  11.  Base  drag  reduction  factor  versus  time  of  flight  for  round  number  4216  fired 
with  propelling  charge  M203E2,  8R,  at  a  quadrant,  elevation  of  499  mils. 
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;  reduction  factor  versus  time  of  flight  for  round  number  4329  fired 
lelline  charee  M203E2,  8R,  at  a  quadrant  elevation  of  748  mils. 
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Figure  13.  Base  drag  reduction  factor  versus  time  of  flight  for  round  number  4219  fired 
with  propelling  charge  M203E2,  8R,  at  a  quadrant  elevation  of  1147  mils. 


r 


Z*0  9*0  S’O  ro  £-0  2*0  1*0  0*0 


s:^napgjao3  S'biq  DinreiiXpoiaY 


21 


Mach  Number 

Figure  14.  Aerodynamic  drag  force  coefficients  for  the  155mm,  DPICM,  M864  base-bum 
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Figure  16.  Factor  fi  versus  Mach  number  for  the  155mm,  DPICM,  M864  projectile. 
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Residuals  of  the  stepwise  multiple  reeression  analysis  versus  Mach  number. 
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List  of  Symbols 


Symbol 

Cdo 

Cdb, 

Cd. 

d 

fo 


h 


/z 

/s 


h 


9 

m 

nir 

P 

Pr 


rt 


rt 

r 

r 


rt 


rt 


Definition 

zero  yaw  drag  force  coefficient 

inert  base-burn  motor,  base  drag  component  of  zero  yaw  drag 
force  coefficient 

radar  determined  drag  force  coefficient 
reference  diameter  of  projectile 

factor  as  a  function  of  quadrant  elevation,  previously  identified 
as  i ,  used  for  matching  experimental  range  firing  data 

factor  used  to  represent  the  drag  reduction  as  a  function  of 
Mach  number 

factor  used  to  represent  the  drag  reduction  as  a  function  of 
time  of  flight 

factor  used  as  a  parameter  for  matching  experimental  range 
firing  data 

factor  used  to  represent  the  drag  reduction  due  to  an  apparent 
effect  of  local  air  pressure 

acceleration  due  to  gravity 

fuzed  projectile  mass  at  time  t 

reference  fuzed  projectile  mass 

air  pressure  at  trajectory  position 

reference  (standard)  air  pressure  (1013.25  mb) 

yaw  of  repose  drag  term  in  the  “Modified  Point  Mass  Trajectory 
Model” 

trajectory  estimated,  slant  range  magnitude 
trajectory  estimated,  slant  range 

HAWK  radar  determined,  rate  of  change  of  slant  range  with  time 

time  derivative  of  the  HAWK  radar  determined,  slant 
range  rate  of  change 

trajectory  estimated,  rate  of  change  of  slant  range  with  time 

trajectory  estimated,  time  derivative  of  the  slant  range  rate  of 
change 
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List  of  Symbols  (Continued) 


Symbol  Definition 

Ut  trajectory  estimated,  velocity  of  the  projectile  with-respect-to  the 

ground-fixed  axes 

Ut  trajectory  estimated,  acceleration  of  the  projectile  with-respect-to 

the  ground-fixed  axes 

Ur  velocity  of  the  projectile  with-respect-to  the  ground-fixed  axis, 

determined  from  HAWK  radar  data  and  estimated  trajectory 

Ur  acceleration  of  the  projectile  with-respect-to  the  ground-fixed  axes, 

determined  from  HAWK  radar  data  and  estimated  trajectory 

speed  of  projectile  with-respect-to  air 

V  velocity  of  the  projectile  with-respect-to  air 

w  velocity  of  the  air  with-respect-to  the  ground  (wind  velocity) 

A  acceleration  due  to  Coriolis  effect 

p  density  (specific  mass)  of  air 
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